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Abstract 



We consider a general non-linear sigma model coupled to Einstein gravity and 
■ show that in spherical symmetry and for a simple realization of self-similarity, 

the spacetime can be completely determined. We also examine some more 
. specific matter models and discuss their relation to critical collapse. 
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INTRODUCTION 



The problem of gravitational collapse has received a great deal of attention in the past 
several years. This is due in large part to Choptuik's results on the collapse of a real scalar 
field coupled to gravity [Q. Choptuik showed that the strong field regime exhibited rather 
striking and unexpected behavior. This included the surprising result that black holes can 
form with arbitrarily small mass and that a scaling relation exists for the black hole mass. 
In addition, the "critical" solution exactly at the threshold of black hole formation possesses 
a discrete self-similarity and serves as an intermediate attractor for near-critical solutions in 
the space of initial data. 

Considerable effort has been made to understand these results from both computational 
and analytic perspectives. The general behavior is now rather well understood and a very 
nice review article summarizes the relevant literature and current state of the art 0. How- 
ever, as pointed out there, progress in understanding critical phenomena from a purely 
analytic perspective has been slower in coming. This is largely due to the well known dif- 
ficulties associated with finding closed form solutions to Einstein's equations in dynamical 
situations. One of the few closed form solutions to the Einstein equations which appears 
somewhat related to critical phenomena is the so-called Roberts solution. It had been writ- 
ten down several times in the literature || Q], but was rediscovered in an attempt to get 
an analytic grasp of the features seen in critical collapse ||. The solution is spherically 
symmetric and continuously self-similar. It also contains timelike singularities. When these 
singularities are cut out and replaced with Minkowski space, the spacetime describes a col- 
lapsing ball of scalar field. In addition, as a parameter is varied, the solution describes 
dispersal of the field to infinity or the formation of an apparent horizon. 

Our attempt here is to extend this solution to a wider class of matter models. Numerical 
work has shown that critical behavior arises in a variety of matter models coupled to gravity. 
Thus one might expect the Roberts solution to be a member of a larger class of these "toy 
models." 

Indeed, we find that the basic spacetime for scalar field collapse as described by the 
Roberts solution turns out to extend virtually unchanged to a large class of massless inter- 
acting scalar fields. It is this result that we present in the next section and follow with some 
specific examples. 



We are interested in considering a generalization of the Einstein-scalar field problem. So 
we turn our attention to Einstein gravity coupled to multiple interacting scalar fields. The 
particular form of the theory we will consider is given by the action 



THE MODEL AND THE SPACETIME SOLUTION 




(1) 



where the Lagrangian density for the matter fields is given by 
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where we have N scalar fields <p A (A = 1 . . . N) and Gab is a set of functions of the scalar 
fields which must be specified to fix the model. In general, the functions making up Gab 
take the form of a Riemannian metric on the iV-dimensional internal space of the scalar 
fields (the target space). Thus the N real scalar fields, 4> A , are coupled to Einstein gravity 
with A a dimensionless coupling constant. In fact, any parameters entering the target space 
metric Gab will also be dimensionless and this in turn will allow solutions to be self-similar. 
The theory defined by Eq. (^) is thus a particular case of more general Einstein-multi-scalar 
theories of gravity such as those discussed in || . 

Varying the action with respect to the metric and the matter fields yields the Einstein 
field equations and the equations of motion for the matter: 

R pv = XG A b(<P C )<P A ,,<P B ,u, (3) 

^V-^.^r^, (4) 

where, as usual, a comma and a semicolon denote partial and covariant differentiation, 
respectively. We also define the Christoffel symbols on the target space 



1 bc - 2 u 



9Gbd dGcD 9Gbc 



(5) 



Because of the contracted Bianchi identities T^ u .\ g uX = 0, Eqs.(|3]) - (H) are not independent. 
In the following, we shall use this fact to solve these equations by properly choosing the 
independent ones. 

Working in spherical symmetry, we can write the metric in double null coordinates 

ds 2 = -2e 2a dudv + r 2 (d6 2 + sin 2 0d<p 2 ) , (6) 

where a and r are functions of the two null coordinates u and v only. The metric is invariant 
under the coordinate transformations u = u(u), v = v(v). In addition to assuming spherical 
symmetry, we are interested in spacetimes that are continuously self-similar (CSS). Stated 
invariantly, we are assuming the existence of a homothetic Killing vector field £ satisfying 



^9ij,u — tiw + C^sM — 2g^, (7) 

where C denotes the Lie derivative. It is then straightforward to show that the homothetic 
vector field is 

9 d 

In addition, it follows that the the metric coefficients take the form 

a = a(z), r = us(z), (9) 

where z = v/u. We also need to know the form of the scalar fields in the presence of the 
continuous self-similarity. The simplest ansatz is to assume 

C^ A = (10) 
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such that the scalar fields will depend only on the coordinate z = v/u as with the metric 
coefficients. This is, in fact, the ansatz we will make here, but it is not necessarily the most 
general. Depending on the symmetries of the N- dimensional target space defined by the 
scalar fields, one possible generalization is to allow the homothetic motion to be accompanied 
by a change in the scalar fields which respects the target space symmetries. For example, 
in the case of a complex scalar field (N = 2 and Gab = Vab, a flat target space), there is 
a U(l) symmetry which could be incorporated into the self-similarity ansatz for the fields 
||. Another possible generalization is to make a kind of hedgehog ansatz for the scalar 
fields before imposing self-similarity. In this case, the (p A fields will have some nontrivial 
angular dependance. An example is the self-similar solution of a spherically symmetric 0(4) 
non- linear sigma model which models the relativistic collapse of a global texture ||. 
For the ansatz given by Eg.(^), the Einstein field equations (|3|) take the form 



2za'(s - zs') + z 2 s 



2 „// 



Xz 2 s 



r< J.A/J.B/ 



'ABS 



As 



2a' V - s" = —G A b$ 



Ai iBi 



2(a' + za") + 2z- 



-,2V 



—XzGab< 

~Z^' 



lAi iBi 



(11) 

(12) 
(13) 

(14) 



where a prime indicates differentiation with respect to the similarity variable z. From 
Eqs.(|i~T|) and ( |i~2"D we find that a = const. By properly rescaling the null coordinates u and 
v, we can set a = 0. Thus, without loss of generality, in the following we shall consider only 
the case where a = 0. Substituting into Eq. (|1"4]) we find that 



(z) = az 2 — z + j3, 



(15) 



where a and (3 are two integration constants. Meanwhile, Eqs. flTTD - ([T3|) reduce to the single 
equation 



G 



ABS 



kAi iBi 



1_ 
4' 



(16) 



where 7 = (1 — Aaf3) /(2\). It is interesting to note that without having specified the target 
space metric Gab-, we have been able to completely determine the metric coefficients. Thus, 
spherical symmetry and our simple ansatz for self-similarity completely fix the spacetime 
geometry to be the same for all of these non-linear sigma models. In particular, the spacetime 
and its global structure is the same as that of the Roberts solution. 



Strictly speaking, this will be true provided 1 — 4a/3 > as in the Roberts solutions. However, 
in the following, we will encounter a situation where this is not necessarily the case. 
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SOLUTIONS OF THE MATTER EQUATIONS 



To complete the solution, we need to specify the non-linear sigma model by choosing a 
form for Gab- In the following we will restrict ourselves to a 2-dimensional target space of 
the form 



Gab 



1 

m 



(17) 



where, for clarity, we define (f) 1 = <fi and </> 2 = ip. 

The matter field equations Eq.(HI) can now be written as 



,2i/ 



s 2 f{4>W 



2 
0. 



■SAM 
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Clearly, the second equation has the first integal, 

dip c 



dz /(0) S 2 ' 

where Co is an integration constant. Inserting the above equation into Eq.( 



d ± = ± 

dz 
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1/2 



(18) 
(19) 

(20) 
we obtain 

(21) 



To solve Eqs. (|20|) and (pT|), we need to specify the function /(</>). In the following, we shall 
consider some particular cases. 

I f{f) = e~+ 

This case corresponds to the non-linear cx-model considered in || for the case k > 0. In 
fact, if we set 



1 1 + ir 



T 



then the Lagrangian density of the complex matter field F takes the form 





VF\ 


2 


(1- 


- K 


F 


\2\2 



16« v 



which corresponds to the one given by Eqs.(Q) and (|T7|) with 



A = -, f{4>) = e~*. 

OK 



(22) 



(23) 



(24) 



in |g with 



As noticed in 0, this also corresponds to the Brans-Dicke (BD) theory of gravity considered 

A = o/ + ! (25) 
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where u is the BD coupling constant. Note the difference between A used here and the one 
used in [|10|1 . 

Substituting Eq.(|24|) into Eqs. (j2"0"|) and ([H]), and then integrating them, we obtain 



21n^ + ln4 



1 + 9' 



c z 



^/7 4 

co 1 + 9 



(26) 



where 



9 = ci 



z — z 4 



z — z^ 



z± 



1 ± VI - 

2^ ' 



(27) 



and Ci, 0o an d V'o are integration constants. Note that in this case the parameters a and /3 
have to be chosen such that af3 < 1/4. Otherwise, no solutions exist. 

II /(0) = sin 2 20 

For this case we have the non-linear a-model considered in || for k < 0. Actually, if we 

set 



tan i 



(sin 2ib + icos2ip) 



(28) 



we find that 





VF 


i 


(1- 


- K 


F 


2)2 



■r-r U' a Q + sin^ 2<pib> a ib c 
\k\ v 



(29) 



which corresponds to Eqs.(@) and (|TT| ) with 



A = — , /(</>) = sin 2 
\k\ 



(30) 



Inserting the above expressions into Eqs.([Hj) and (21) , we find the following solution 

1/2 



cos 2<p 



2Ac 2 



l-Aa(3 



COS£f, 



1 



ib = +— tan 1 

^ 2 



tan (yf 

Co 



(31) 



but now 



(jf = ±a/|k| In 



,2 — Z 4 



z — z_ 



Cl, 



(32) 



where zj- are given by Eq . (^7|) , and ci and are other integration constants. Similar to 
the last case, now the condition a(3 < 1/4 has to hold. 
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Note that when k = 0, the self-similar solution of the non-linear cr-model of Eq.(p3|) or 
Eq.(H) with the ansatz (jT^) is the same as that of Roberts Q. As a matter of fact, in the 
latter case we have f(<f>) = 1. Then, from Eqs. ( p0|) and (|2T| ) we find that ^ = A/>, where A 
is a constant. 



ni f{4>) 



A = ~ (3 + 2uj) 



Finally, we consider the case corresponding to the BD theory considered in JT0[ . As noted 
previously, when uj > —3/2, it is equivalent to the non-linear a-models with k > 0. However, 
when uj < —3/2, such a correspondence no longer exists. Therefore, in the following we shall 
consider only the case where uj < —3/2. Before proceeding and by way of justification, we 
note that the case where uj < is also worthy of study for more than its intrinsic theoretical 
interest. In particular, for sufficiently negative uj, the theory passes all the experimental and 
observational constraints as does the case of positive uj ||1 1|| . However, when uj < —3/2, the 
static spherical spacetimes are quite different from the ones with uj > —3/2. In particular, 
in the former case "cold" spherical black holes exist [12], while in the latter case, the only 



static spherical black hole is the Schwarzschild black hole [fL3|| . 

Inserting the above expressions for A and f(<j)) into Eqs. 02 
solutions exist for both a [3 < 1/4 and a [3 > 1/4. 

Ilia 1 - Aa(3 > 
In this case, we have 



and (Elf), we find that 



In 



cos 2 g 



+ 



if) = Cq|3 + 2uj\ tang + xjj 



(33) 



where 



and z± are given by Eq.(|2 

Illb 1 - Aa(3 < 
Now we have 



In 



z — Za 



~ (8A 2 cg)V2 
and e = sign{ot). 



-In 



sinh 



+ 



where 
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± 



2cI(2uj + 3) 


1/2 

tan -1 


laP- 1 





Cl, 



(34) 



: - a ^ ^-cotanh ( — + ib Q , 

c 2a; + 3 V 2 J 



2az - 1 



(35) 



(4a/3 -1)V2 



+ Ci. 



(36) 



For r to be real, we need to impose the condition a > and (3 > in this case. 

Clearly, by giving different /(</>), we can obtain various solutions for the matter fields 4> 
and ip. However, these solutions are not all independent. For some particular choices, they 
may correspond to the same matter fields. This follows, of course, from the diffeomorphism 
invariance of the target space manifold. For example, the choice /(</>) = sinh 2 2<p gives the 
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same solution as that given in Case I. This is closely related to the fact that the matter 
fields described by the complex function F has a SL(2, R) symmetry for the case k > ||. 
But for k < 0, no such symmetry exists. We would thus expect that the solution given in 
Case II is unique. 

Since the spacetime geometry is independent of the particular matter model, the global 
structure of the spacetimes should be the same as that of the Roberts solution ||, provided 
that 1 — 4a (3 > 0. Of course, the Ricci scalar and the mass will also be unchanged 

7/7; 7/7; 

R=^(l-4a(3) M = ~(l-4ap). (37) 

Again, one can remove the timelike singularities which are present and replace them with 
Minkowski space, matching at the advanced time v — 0. In this way, we imagine the scalar 
fields to be turned on at v — 0. Thus, the solution with /3 — 1 resemble critical collapse 
in the sense that for a < black holes are formed (supercritical regime), while for the case 
< a < 1/4, no black holes are formed and the fields disperse to infinity (subcritical regime) 
||. Note, however, that because of the matching, the solutions are C° in contrast to the 
smooth, attracting critical solutions which have been found in numerical calculations. 

When 1 — 4a/3 < 0, which is the case Illb, there is no correspondence to the Roberts 
solution. From Eq . fl37|) we can see that to have the mass be non-negative, the solution has 
to be restricted to the regions u > 0, v > and u < 0, v < 0, which will be referred, 
respectively, as regions / and //. In region II the solution can be considered as representing 
the gravitational collapse of two scalar fields and a point-like singularity is finally formed 
at u — = v. This singularity is direction-dependent. When the singularity is approached 
along different directions, the singular behavior will be different. In particular, if one ap- 
proaches the singularity along the null hypersurface u = or the one v = 0, everything seems 
regular. On the hypersurfaces u — 0, v < and v — 0, u < 0, the two scalar fields carry null 
data. Using the arguments for the case 1 — 4a(3 > 0, we can replace the solution in the 
regions u > 0, v < and u < 0, v > by two flat regions. Once this is done, the spacetime 
needs to be further extended beyond the hypersurfaces v — 0, u > and u = 0, v > 0, 
which are Cauchy horizons. The extension beyond them are not unique. For example, one 
can match the two flat regions to the solution given by Eq.(|5o]) for u > 0, v > 0, or make 
an analytic extension and replace region I by another flat region. However, in either case 
the singularity at u = = v is naked. Figure 1 gives a conformal diagram of the spacetime 
summarizing the above results. 

Finally, we would like to note that the conclusion that the spacetime geometry is in- 
dependent of the choice of the particular scalar field model is based on the fact that the 
spacetime realizes a fairly simple assumption of self-similarity. Otherwise, it is easy to show 
that this conclusion does not hold. Therefore, for the general perturbations (which do not 
have such a symmetry), it would be expected that they will depend on the choice of non- 
linear sigma model. In particular, there are no a priori reasons to expect that the spectrum 
of the unstable modes of the perturbations of the solution a = and (3 = 1 for a general 
choice of Gab is the same as that for the particular Roberts solution as was considered in 



14]. This point deserves further investigation and will be considered in future work. 
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FIGURES 

FIG. 1. The conformal diagram for the spacetime described by case Illb. The singularity is 
at (u, v) = (0, 0) and the Cauchy horizons are at u = 0, v > and v = 0, u > 0. Region / can be 
replaced either (analytically) by flat space or (continuously) by the solution given by Eq. (|35|) . For 
more details, see the text. 
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